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The  hazard  posed  by  trailing  vortex  wakes  of  Large  aircraft 
has  been  of  increasing  concern  due  to  its  effects  on  flight  safety 
and  because  of  the  economic  impact  of  the  required  increase  in 
aircraft  separation  during  takeoff  and  landing.  Investigations  of 
vortex  wake  behavior,  both  experimental  and  analytical,  are  being 
employed  to  set  guidelines  for  avoidance  procedures  and  to  deter- 
mine methods  of  alleviation.  Purely  analytical  methods  have  had 
difficulty  achieving  agreement  with  measured  vortex  behavior  due 
to  the  complex  flow  processes,  such  as  turbulence,  occurring  with- 
in and  about  the  vortex.  Since  vortex  descent  within  the  atmos- 
phere introduces  the  added  complexity  of  density  stratification, 
numerical  methods  must  be  used  to  obtain  any  hope  of  including 
the  important  physics  of  the  problem. 

This  investigation  of  vortex  descent  has  been  undertaken  to 
clarify  the  results  of  various  analytical  studies  which  have  been 
made  by  numerous  investigators  and  which  have  led  to  contradictory 
predictions  of  vortex-pair  behavior  in  a stratified  atmosphere. 

The  calculations  which  are  reported  here  are  based  on  the  solution 
of  the  mean  momentum  and  Reynolds  stress  equations  of  turbulence 
based  on  a second-order  closure  model.  The  computer  code  which  pro- 
vodes  this  solution  has  been  given  the  designation  WAKE  . To  pro- 
vide verification  of  this  code  for  the  purpose  of  this  study,  vor- 
tex rings  descending  in  a stratified  fluid  are  measured  experimen- 
tally, and  comparisons  are  made  between  results  of  the  axisymmetric 
version  of  the  code  and  data  for  the  conditions  of  the  tests. 

The  flight  data  selected  for  comparison  and  analysis  are  the 
measured  vortex  descent  generated  by  a Boeing  747  in  landing  con- 
figuration, but  in  level  flight.  The  code  is  run  using  an  initial 
vortex  rotational  velocity  profile  which  is  close  to  the  actual 
data.  Because  the  initial  vortex  is  so  strong  and  tight. 


a new  method  has  been  devised  for  calculating  the  vortex  behavior 
based  on  perturbation  about  the  initial  condition. 

The  results  to  be  presented  here  are  built  upon  past  efforts 
performed  for  several  agencies  which  deserve  mention  and  acknowledg- 
ment. Development  of  the  turbulent  model  has  been  supported  in  large 
part  by  the  Environmental  Protection  Agency,  under  Contracts  68-02- 
0014  and  68-02-1310.  The  WAKE  computer  code  has  been  developed  in 
its  two-dimensional  form  through  funding  from  these  contracts  and 
through  support  by  the  Defense  Advanced  Research  Projects  Agency, 
as  monitored  by  the  Office  of  Naval  Research  \inder  Contract  N00014- 
72-C-0413.  The  axisymmetric  version  has  been  developed  primarily 
through  efforts  performed  for  NASA  Ames  on  Contract  NAS2-9037  and 
for  the  Nuclear  Regulatory  Agency  xander  Contract  AT(49-24)-0218. 

The  influence  of  stratification  on  vortex-pair  behavior  has  been 
studied  with  the  code  under  contract  to  NASA  Langley  under  Contract 
NASl-14707. 

The  tank  and  much  of  the  apparatus  required  to  perform  the 
vortex  ring  experiment  have  been  fabricated  under  funding  from 
the  John  Hopkins  Applied  Physics  Laboratory,  Contract  APL-600646. 

The  contract  monitor  for  this  investigation  is  James  N.  Hallock. 
The  authors  gratefully  acknowledge  his  help  and  that  of  David  C. 
Burnham  for  providing  sources  of  flight  data  during  the  data 
selection  process  of  this  effort. 
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V vortex-ring  or-pair  descent  velocity 

w core  swirl  velocity 

z maximum  descent  distance  of  vortex  ring 

max 

r vortex-ring  or-pair  circulation 

e atmospheric  dissipation  rate 

t vorticity 

n radial  coordinate  measured  form  center  of  vortex 

ring  core 

0 mean  temperature 

0 fluctuating  temperature 

A turbulent  macroscale  parameter 

X turbulent  microscale  parameter 

V kinematic  viscosity 

p density 

stream  function. 


Superscripts 

(~)  quantity  normalized  by  r. 


Subscripts 

avg 

c 

max 

o 

P 

Special 

<> 


average 
countersign 
maximum  value 
initial  value 

refers  to  principal  vortex, 
ensemble  averaged. 
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Coordinate  System 


Cartesian  Coordinates 


u.v.w  velocities  parallel  to  x,y.  and  r coordinates, 

respectively . 

x.y.z  coordinates  parallel  to  vortex  pair  cores, 

horizontal,  and  vertical,  respectively. 

Circular  Cylindrical  Coordinates 

r,0,2  coordinates  In  radial  (horizontal),  azimuthal, 

and  vertical  directions,  respectively 

u,v,w  velocities  In  the  r,0,  and  z directions,  respectively. 
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1.  INTRODUCTION 


1.1  VORTICES  IN  AN  INVISCID  FLUID 

The  hazard  created  by  aircraft  trailing  vortices  has 
generated  renewed  interest  and  research  into  the  physics  of 
vortex  flows.  In  an  ideal  fluid  without  external  influences, 
vortex  pairs  and  vortex  rings  will  persist  indefinitely.  For 
a ring  with  uniform  concentrated  vorticity,  the  speed  of  trans- 
lation has  been  known  for  over  a century.  Kelvin,  in  1867, 
gave  the  translation  velocity  for  a vortex  ring  in  an  ideal 
fluid 


[■■(if)-i] 


(1.1) 


where  R is  the  ring  radius,  and  r^  is  the  core  radius  within 
which  the  vorticity  is  uniform.  The  self-induced  velocity  of 
translation  of  vortex  rings  with  arbitrary  distribution  of  vor- 
ticity in  the  core  was  determined  independently  by  three  inves- 
tigators only  recently:  Bliss,  Ref.  1,  Saffman,  Ref.  2,  and 
Fraenkel,  Ref.  3.  The  velocity  of  translation  was  found  to  be 


is  the  core  radius  within 


(6  ^ ■ I • 


(1.2) 


where  A is  defined  by 

. tim 
A = n/r  - 


P (")  J ■ ■ 


(1.3) 


ana  v is  the  swirl  velocity  non-dimensionalized  by  r/4TTR, 
and  n is  a radius  measured  from  the  center  of  the  core  ( see 
Fig.  1-la).  For  the  case  of  a core  with  uniform  vorticity, 

A = ^ for  r^/R<<l,  recovering  Eq.  (1.1). 
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For  a vortex  pair  with  concentrated  vorticity, the  speed 
of  translation  is  not  critically  dependent  upon  the  distribu- 
tion of  vorticity  within  the  core.  For  a vortex  pair  with  con- 
centrated vorticity  in  the  cores  with  spread  r^  , the  velocity 
for  core  separation  b is 

The  coordinates  for  a vortex  pair  are  shown  Figure  (1-lb) . 


1.2  VORTICES  IN  REAL  FLUIDS 
1.2.1.  Vortex  Pairs 

In  the  presence  of  buoyancy,  shear,  or  three-dimensional  in- 
stabilities, the  life  time  of  vortex  pairs  or  vortex  rings  will 
be  limited  by  one  or  more  of  these  effects,  and  the  velocity  of 
translation  will  no  longer  obey  Eqs . (1.2)  or  (1.4). 

In  a real  fluid,  the  action  of  viscosity  and  turbulence 
will  restrict  the  distance  vortices  attained  during  descent  to 
finite  values.  The  study  of  aircraft  trailing  vortices  must, 
therefore,  include  the  influence  of  all  of  these  effects,  which 
occur,  or  may  occur, during  the  descent  of  vortices  in  the  real 
atmosphere . 

The  role  of  atmospheric  stratification  and  buoyancy  in  in- 
fluencing the  behavior  of  trailing  aircraft  vortices  has  been 
the  focus  of  considerable  attention  in  the  recent  past.  Table 
1-1 , adapted  from  Widnall,  Ref.  4,  presents  a summarized  compari- 
son of  the  predicted  behavior  of  vortex  pairs  in  a stably  stra- 
tified environment.  Refs.  5-12.  Appended  to  Widnall 's  compari- 
son is  the  recent  work  of  Bilanin  et  al..  Ref.  13.  There  is 
clearly  no  consensus  regarding  vortex  spacing,  descent  velocity, 
and  the  effect  of  buoyancy-generated  vorticity.  All  investiga- 
tors, except  the  last,  disregard  the  effects  of  viscosity  and  tur- 


3 


TABLE  1-1.  COMPARISON  OF  THEORETICAL  MODELS  FOR  THE  DESCENT 
OF  A VORTEX  WAKE  IN  A STABLY  STRATIFIED  ATMOSPHERE 


Author  ( Ref. ) 


V ortex 
S pactn 


Descent 
S peed 


Buoyancy- 
G enerated 
V orticit 


Costen  (5)* 

1 

1 

Weak  stability 

Strong  stability 

Decreases 

Increases 

Increases 

Stops 

No  effect 

Kuhn  & Nielsen  (6) 

Decreases 

Increases 

Partly  entrained 

Saffman  (7) 

Increases 

Stops 

Entrained 

Scorer  & Davenport  (8) 

Decreases 

Increases 

Detrained 

Tombach  (9) 

Strong  stability 

Weak  stability 

Increases 

Decreases 

Stops 

Increases 

Entrained 

Entrained 

Tulin  & Shwartz  (10) 

Increases 

Stops 

Entrained  and 
annihilated 

Crow  (11) 

Decreases 

Increases 

Created  in  wake 
and  interface 

Narain  & Uberoi  (12) 

Decreases 

Decreases 

Entrained 

Bilanin  et  al.  (13) 

Approx . 
constant 

Decreases 

Complex 

interaction 

*Costen's  model  is  a buoyant-core  model;  all  other 
oval  models 

s are  buoyant- 

bulent  transport  in  their  analytic  model  although  all  appear  to 
agree  on  the  importance  of  these  effects.  There  is  no  dispute 
regarding  the  importance  of  viscosity  and  turbulent  transport 
in  controlling  the  behavior  of  descending  vortex  pairs,  but  only 
in  how  they  enter  into  the  physics  of  the  flow  and  if  any 
proposed  predictive  model  accounts  for  the  physics  correctly. 

1.2.2  Vortex  Rings 

In  recent  years,  there  has  been  considerable  effort  applied 
to  the  understanding  of  vortex  rings  although  much  of  it  has 
not  been  concerned  with  vortex-ring  behavior  in  a stratified 
fluid.  Turner,  Ref.  14,  presented  a theoretical  study  of  the 
comparison  between  buoyant  vortex  rings  and  vortex  pairs. 

Linden,  Ref.  15,  conducted  experiments  and  proposed  a theoretical 
model  for  the  interaction  of  a vortex  ring  with  a sharp  density 
interface.  Linden's  experiments  were  conducted  in  the  low  Froude 
number  regime  and  instead  of  penetrating  completely  into  the 
high-density  region,  the  rings  rebounded  at  the  density  inter- 
face . 

Recently,  Maxworthy, Ref . 16,  discussed  results  of  stratified 
vortex  ring  experiments . The  experiments  were  carried  out  in  a 
linearly  stratified  environment,  and  the  vortex-ring  radii  were 
observed  to  decrease  with  increased  penetration  into  the  strati- 
fication. Maxworthy  discussed  the  experimental  behavior  within 
the  context  of  a broader  model  of  vortex  ring  behavior  developed 
in  Refs.  16-18.  However,  as  noted  by  Maxworthy,  the  complex 
interactions  occurring  among  the  ring  dynamics,  buoyancy,  and 
turbulent  transport  cannot  be  fully  accounted  for  by  a simple 
model.  It,  therefore,  seems  clear  that  we  must  resort  to  numer- 
ical calculations  which  contain  a more  complete  description  of 
these  interactions. 
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1.3  V/\LIDATION  OF  THE  PREDICTIVE  MODEL 


The  test  of  any  proposed  model  is  whether  it  predicts  the 
observed  behavior  accurately.  There  are  full-scale  data  available, 
but  little  data  exist  which  have  either  been  taken  under  controlled 
conditions,  contain  the  relevant  measurements  sufficiently  to 
validate  predictive  models,  or  are  without  ground  interference 
effects.  Although  laboratory  experiments  seem  to  offer  the 
logical  alternative  to  obtain  experimental  verification  of  theory, 
it  is  difficult  if  not  impossible  to  study  in  the  laboratory  a 
vortex  pair  free  of  end  effects  (e.g.,  wall  scrubbing)  which  in- 
fluence their  behavior. 

The  study  of  vortex  dynamics  in  the  laboratory  is  most  easily 
accomplished  by  observing  the  behavior  of  vortex  rings.  This 
approach,  therefore,  offers  an  alternative  in  predictive  model 
evaluations  and  in  increasing  our  understanding  of  vortex  behavior. 

The  following  three  sections  of  this  report  contain  the  re- 
sults of  such  an  experiment,  and  present  a comparison  between 
data  and  the  results  of  calculations  made  with  an  axlsymraetrlc 
stratified  code.  This  code  incorporates  the  effects  of  turbulence 
based  upon  second-order  closure  of  the  Reynolds  stress  equations, 
using  an  invariant  model  of  turbulent  transport  (summarized  by 
Donaldson,  Ref.  19;  Lewellen,  Ref.  20;  and  Level len,  Teske,  and 
Donaldson,  Ref.  21).  The  model  has  been  incorporated  into  both 
two-dimensional  and  axisymmetric  computer  codes  which  have  been 
used  with  considerable  success  in  the  study  of  a variety  of 
stratified  wake,  aircraft  trailing  vortex,  and  meteorological 
phenomena,  see  Refs.  13  and  22-25.  These  codes  are  described  more 
fully  later  in  this  report. 

The  vortex  t ing  experiment  and  apparatus  are  described  in 
Section  2.  In  Section  3,  the  axisymmetric  WAKE  code  is  described, 
and  the  vortex  ring  experiment  initial  conditions  are  discussed. 
Section  4 presents  comparisons  between  the  vortex  ring  data  and 
numerical  calculations.  Tl\c  two-dimensional  WAKE  code  is  briefly 
described  in  Section  5,  the  vortex  pair  flight  data  selection  is 
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described  and  comparisons  between  numerical  and  flight  test  results 
are  made.  Conclusions  are  drawn  and  discussed  In  Section  6. 


2.  DESCRIPTION  C?  THE  VORTEX  RING  EXPERIMENT 


2.1  EXPERIMENTAL  APPARATUS 

The  apparatus  used  to  perform  the  vortex-ring  experiments  in 
a stratified  fluid  consisted  of  a transparent  test  tank,  a nozzle- 
type  vortex-ring  generator,  and  the  necessary  mixture  tanks,  pumps, 
and  piping  required  to  produce  a density  stratified  fluid.  A 
sodium  chloride/water  brine  solution  constituted  the  test  medium. 

2.1.1  Test  Tank 

The  tank  was  constructed  to  permit  observation  of  the  ring 
formation  and  descent  and  the  use  of  quantitative  flow  visuali- 
zation techniques . A schematic  of  the  tank  is  shown  in  Figure 
2-1.  The  tank  had  inside  dimensions  of  142  cm  in  depth  and  was 
square  with  sides  of  61  cm.  Two  opposite  sides  and  the  bottom 
of  the  tank  were  made  of  plywood  and  covered  by  a white  vinyl 
liner.  The  remaining  sides  were  of  1.9  cm  acrylic  sheet.  A 
combination  fill-drain  fitting  was  installed  at  the  center  of  the 
bottom.  To  reduce  concentrated  disturbances  caused  by  the  filling 
process,  a screen  diffuser  was  mounted  2 cm  above  the  bottom. 

Filled  to  within  7 cm  of  the  top,  the  tank  held  500  liters,  or 
about  410  kg  of  fresh  water. 

2.1.2  Water  Feed  System 

The  system  used  to  produce  the  stratified  brine  solution  is 
shown  schematically  in  Figure  2-2.  As  described  by  Oster,  Ref. 26, 
a linear  gradient  is  produced  by  setting  the  rate  of  flow  of  the 
brine  supply  at  one-half  of  the  rate  of  flow  of  the  mixture  out 
of  the  mixing  tank,  which  initially  holds  fresh  water.  The  initi- 
al voltame  of  fresh  water  is  specified  by  the  requirement  that  the 
mixture  tank  run  dry  at  the  desired  fill  level  of  the  test  tank. 

The  mixture  tank  is  fitted  with  paddles  to  insure  uniform  mixture 
concentration.  Volume  flow  rate  of  brine  and  mixture  are  measured 
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142  cm 


Figure  2-1.  Vortex— ring  experimental  tank. 
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by  a pair  of  calibrated  flow  meters. 

2.1.3  Vortex  Ring  Generator 

The  device  used  to  produce  vortex  rings  is  shown  in  cross 
section  in  Figure  2-3.  Three  adjustments  were  employed  to  produce 
a wide  variety  of  ring  starting  conditions.  The  length  of  the 
stroke  was  controlled  by  the  positions  of  the  travel  stop  nut 
and  the  trigger  setting  (pull  back)  nuts.  The  speed  of  the 
stroke  was  controlled  by  the  main  spring  tension  and  also  by  the 
trigger  setting.  Once  a particular  set  of  adjustments  was  made, 
the  ring  generator  produced  Identical  rings  repeatedly.  The  acrylic 
nozzle  had  a diameter  of  7.62  cm.  Dye  was  injected  prior  to  each 
run  for  flow  visualization  through  a manifold  which  had  24  equally 
spaced  holes. 

By  tripping  the  trigger,  a cylindrical  slug  of  fluid  is 
ejected  from  the  end  of  the  nozzle.  As  the  fluid  emerges  from 
the  nozzle,  there  are  two  effects  which  influence  the  ultimate 
characteristics  of  the  ring.  First,  there  is  the  weak  ring  of 
vorticity  of  opposite  sign  which  forms  due  to  entrainment  of  the 
fluid  external  to  the  nozzle  (see  Maxworthy,  Ref  16).  Second, 
there  is  a "stopping"  effect  as  the  diaphragm  of  the  ring  maker 
comes  to  rest.  It  was  noted  during  the  test  that  there  was  a 
slight  rebound  which  pulled  a small  portion  of  the  ejected  slug 
of  fluid  back  into  the  nozzle.  Neither  effect  was  a serious 
problem,  since  rings  were  produced  which  were  free 
of  such  "starting"  effects  within  about  one-half  to  one  nozzle 
diameter  away  from  the  exit. 


2.2  FLOW  VISUALIZATION  TECHNIQUE 

Two  methods  of  observing  and  recording  the  vortex-ring  be- 
havior were  used.  The  first  made  use  of  the  injected  dye  as  a 
streamline  tracer.  In  this  case,  a sheet  of  ground  glass  with 
a reference  grid  drawn  on  it  was  mounted  outside  the  rear  wall 
of  the  tank  to  serve  as  a diffuser.  Photoflood  lamps  were  direct- 
ed at  a white  reflecting  panel  behind  the  ground  glass  so  as  to 
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Figure  2-3.  Ring  generator. 


aclilove  uniform  Illumination.  With  this  illumination,  the  ring 
atroamlines  couhl  bo  visualized  against  the  uniform  background. 
Since  there  was  considerable  scattering  of  light  within  the  tank, 
all  sides  of  the  ring  were  also  illuminated,  creating  an  enhanced 
tltree-dimensional  effect.  This  technique  was  particularly  useful 
for  observing  the  details  of  the  ring  structure,  since  the  dye 
tended  to  concentrate  in  certain  regions  more  than  others. 

A second  teclinique  which  was  assessed,  but  which  was  not  used 
during  the  present  test  data  reduction,  was  the  shadowgraph  tech- 
nique. In  this  case,  the  ground  glass  was  moved  to  the  front  wall 
of  the  tank  and  tlie  collimated  beam  of  light  projected  from  the 
rear.  Indications  are  that  this  technique  would  be  quite  useful 
in  studying  the  ring  flow  field  as  it  descends  in  the  tank.  Of 
particular  interest  is  visualization  of  the  turbulent  structure. 

2 . 3 PllOTOGRiM’llY 

(R) 

The  behavior  of  the  rings  was  recorded  by  means  of  a Bolex^ 
lO-imn  motion  picture  camera.  Sequences  were  taken  in  slow  motion 
i04  friimes  per  second)  . 

All  measurements  were  scaled  using  the  nozzle  orifice  dia- 
meter as  the  reference  length.  The  amount  of  refraction  produced 
by  the  vertical  density  gradient  was  found  to  be  negligible  for 
all  conditions  except  during  passage  through  the  sharp  density 
dlscontlnvilty . 

2.4  CALIBRvVnON 

To  determine  the  range  of  performance  of  the  ring  generator, 
a series  of  calibration  runs  were  carried  out  with  the  tank  full 
of  fresh  water.  For  this  purpi>se , t lie  main  spring  tension  was 
set  at  four  different  values.  For  each  of  these  values,  the 
trigger  adjustment  (pull  back)  was  also  set  at  four  values.  Thus, 
a basic  set  of  In  conditions  was  established.  For  each  condition, 
five  trials  were  made  in  which  the  ring  descent  was  timed  from  a 
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few  cm  below  the  nozzle  to  within  a few  cm  above  the  tank  bottom. 
The  total  ring  travel  during  timing  was  91  cm.  A plot  of  the 
calibrations  in  terms  of  average  descent  velocity  versus  trigger 
setting  for  each  main  spring  tension  is  given  in  Figure  2>4.  It 
is  seen  that  the  velocity  is  approximately  linear  with  spring  ten- 
sion and  that  a range  of  velocities  from  about  6 to  37  cm/sec  could 
be  achieved. 

2.5  STRATIFICATION  MEASUREMENT 

For  the  linear  gradient  tests,  a hydrometer  consisting  of  a 
set  of  balls  of  known  specific  gravity,  in  an  open  triangular 
cage,  was  used  inside  the  tank  to  assess  the  accuracy  of  the  den- 
sity gradient  produced.  As  shown  in  Figure  2-5,  the  stratifica- 
tion achieved  was  very  nearly  linear. 


Average  descent  velocity,  cm/sec 
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Figure  2-5.  Measured  stratified  density  (o  Measurement, 
linear  gradient) . 
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3.  VORTEX-RING  NUMERICAL  CALCULATION  TECHNIQUE 


3.1  DESCRIPTION  OF  THE  AXISYMMETRIC  WAKE  CODE 

The  equations  governing  the  mean  variables  and  the  turbulent 
second-order  correlations,  in  modeled  form,  and  based  on  the 
Boussinesq  approximation,  are  given  in  Appendix  A (in  an  axisym- 
metric  coordinate  system) . A derivation  of  the  Reynolds  equations 
for  the  second-order  turbulent  correlations  and  a discussion  of 
the  model  development  and  constant  evaluation  may  be  found  in 
Ref.  20.  The  equations  are  finite  differenced  on  a nonuniformly 
spaced  grid  and,  with  the  exception  of  the  equation  governing  the 
mean  density,  solved  using  a modified  alternating-direction- impli- 
cit scheme.  The  mean  density  is  solved  for  using  a flux-corrected 
transport  technique  developed  by  Boris  and  Book,  Ref.  26.  This  is 
done  to  improve  resolution  of  the  sharp  density  gradient  between 
the  descending  oval  and  the  ambient  fluid.  This  gradient  in- 
creases with  time  as  the  ring  penetrates  progressively  farther  into 
the  stratified  environment. 

The  Poisson  equation  governing  the  stream  function  from  which 
the  velocities  are  obtained  is  evaluated  using  a fast  elliptic 
solver  developed  by  Swartztrauber  and  Sweet,  Ref.  28.  A J6 
representation  of  the  convective  operator  (Arakawa,  Ref.  29),  in 
which  the  conservative  and  nonconservative  operators  are  averaged, 
is  used  to  improve  simulation  of  the  convective  nature  of  the 
vortex-flow  problem. 

To  keep  the  vortex  ring  in the  computational  domain,  an  up- 
wash  velocity  is  added.  Tracking  the  centroid  of  the  square  of 
the  vorticity  is  a technique  which  centers  the  ring  in  the  domain 
until  the  very  final  stages  of  the  total  descent. 
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The  computations  were  initialized  by  positioning  a 
Gaussian  of  vorticity  of  the  form  , 


in  the  azimuthal  plane  at  a point  r /R  and  z /R  , which  nor- 

mally  would  be  at  1.0  and  0.0,  respectively.  However,  i’q/Rq  “ 0.95 

was  used  here  for  reasons  which  will  be  described  later.  The 

coordinate  n is  measured  radially  outward  from  the  vortex  core 

center,  and  r^  is  the  vortex  spread.  An  isotropic  Gaussian  spot 

of  turbulence  was  positioned  in  the  center  of  the  vortex  core 
2 2 

with  strength  /V^  = 0.01  such  that 

<uu>  = <vv>  = <ww>  = (q^^/3)exp  [-(n/RQ)^/(r(,/Rg)^j  ■ (3.2) 

The  initial  conditions  are  shown  schematically  in  Figure  3-1. 

At  the  edges  of  the  computational  domain,  except  at  r • 0 , 
the  stream  function  is  evaluated  by  a moment  expansion  of  the 
Biot-Savart  law  as  described  in  Appendix  B.  At  the  inflow 
boundary  at  z_  , the  density  is  convected  in  at  the  appropriate 
value  for  the  type  of  stratification  considered  and  the  descent 
which  has  occurred.  For  a linear  stratification,  the  density  is 
non -dimens ionalized  by  the  product  of  the  background  gradient 
and  R^  , so  that  3p/az  » -1  in  non -dimensional  form.  For  the 
case  of  the  sharp  density  discontinuity,  the  total  non-dimensional 
density  difference  across  the  discontinuity  is  taken  as  unity, 
with  a parameter  used  to  specify  the  vertical  spread  of  the 

density  interface.  The  full  expression  used  for  the  nondimen- 
sional  density  is 


P - -0.5  i|J-  erfdzI/Cg)  (3.3) 

It  was  found  that  * 0.6  provided  a density  profile  which 

was  satisfactory  in  simulating  the  actual  thickness  of  the  density 
interface  in  the  experiment. 
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At  the  outflow  boundary  at  , an  advection  condition  was 

used  to  sweep  the  flow  through  the  boundary.  The  turbulence  and 
vorticity  inflow  values  were  set  equal  to  zero.  Because  of  the 
relatively  high  upwash  velocity, no  wave  interaction  problems 
were  encountered  at  the  r_^  boundary  since  they  were  advected 
out  of  the  domain  before  they  could  reflect  back  into  the  compu- 
tation. 


3.2  INITIAL  VORTICITY  PROFILE 


Because  it  was  not  possible  to  resolve  the  swirl  velocity 

in  and  around  the  core,  another  method  was  devised  to  obtain  the 

appropriate  initial  maximum  vorticity  t and  spread  r 

°max  ^ 

This  was  done  by  utilizing  the  measurement  of  the  oval  semi- 
axis a as  follows.  The  descent  velocity  V of  a vortex  ring 
with  arbitrary  distribution  of  vorticity  in  the  core  is  given  by 
Eq.  (1.2).  For  a Gaussian  distribution  of  vorticity, 

A = -0.058  (Saffman,  Ref  2).  Now  , 


o o 
max 


TR. 


o 


TTr 

c o 


(3.4) 


for  a Gaussian  distribution  so  that  the  descent  velocity  non- 

dtmensionalized  by  V is 

■'  o 


At  tV^/R^  = 0.,  V/V^  = 1.,  R/R^  = 1.,  so  that  if 

r /Rq  can  be  defined,  we  can  obtain  ;; 

max 


(3.5) 


The  semi-axes  a and  b are  defined  by  the  separation 
streamline  ij^(a,0)  = 0 and  tp(0,b)  = 0,  respectively,  of 

a vortex  ring  in  steady  coordinates.  For  a circular  ring  of 
concentrated  vorticity  , 
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.(r.z)  - [(I  - k)  K(k)  -(l)E(k^-  X , 0.6) 

where  K(k)  and  E(k)  are  the  complete  elliptic  integrals  of  the 
first  and  second  kinds,  and  where  V is  given  by  Eq . (1.2) 
with  A • 1 *♦.  k is  defined  by 


Evaluation  of  Eq.  (,3.5)  for  - 0 at  (a.O)  and  (O.b) 
as  a function  of  r^  yields  a and  b as  shown  in  Figure  3-2. 

Both  are  rather  weak  functions  of  r /R  ; , evaluated  using 

max 

Eq.  (3.5)  for  V/V^  - 1 is , however,  a strong  function  of 
r^/R^  for  r^/R„  < 0.3,  Figure  3-3.  Also  shown  is  the  value  of 

Rptp  /V^  vs.  for  a concentrated  core  of  constant  vorti- 

city . 

Using  the  experimentally  determined  values  of  a/R  and 
b/R^  as  a guide  provided  the  means  of  determining  the  core 
spread  ■ This,  with  the  requirement  that  V/V^  = 1.0 

initially,  gave  the  maximum  vorticity. 

As  noted  earlier,  the  vortex  rings  were  initialized  at 
r/R^  “0.95  , “ 0 . This  was  done  to  provide  closer  agree- 

ment with  the  experimentally  observed  initial  radii  of  most  of 
the  runs;  there  was  apparently  still  some  vorticity  rollup 
occurring  at  the  selected  initial  stations. 

The  second-order  turbulent  closure  model  contained  a dynamic 

equation  for  the  scale  A . The  calculations  presented  here  were 

made  for  a constant  scale  length,  and  therefore,  the  scale  equation 

was  not  employed.  A parametric  study  was  carried  out  for  a 

range  of  values  of  A/R  between  0.03  and  0.3.  It  was  found 

o 

that  A/R^“  0.03  gave  closest  agreement  with  data  and  this  value 
for  the  scale  was  used  throughout  the  calculations  to  be  present- 
ed here. 


23 


a/R,  b/R 


4.  COMPARISON  OF  NUMERICAL  RESULTS  WITH  VORTEX-RING  EXPERIMENTAL 

DATA 


4.1  DESCRIPTION  OF  DATA 

Those  experimental  cases  selected  for  analysis  are  listed 
in  Table  4-1.  Also  shovm  are  the  core  spread  and  maximum  vorti- 
city  used  in  initializing  the  WAKE  code  calculations.  Case  1 
was  anon-stratified  baseline  rtin;  cases  2 and  3 were  rings  of 
fresh  water  of  specific  gravity  1.0  injected  into  a brine 
mixture  of  specific  gravity  1.19.  The  density  discontinuity 
was  located  approximately  three  ring  diameters  below  the  nozzle 
exit.  Cases  4 through  6 were  conducted  in  a linearly  stratified 
density  fluid.  The  nozzle  was  immersed  in  fresh  water,  and  the 
linear  stratification  started  at  approximately  two  nozzle  diame- 
ters below  the  nozzle  exit.  This  was  done  to  insure  that  the 
rings  consisted  of  a known  density  since  some  mixing  was  inevi- 
table during  and  shortly  after  each  run.  The  location  of  the 
hydrometer  balls  was  noted  before,  during,  and  after  the  test, 
and  showed  that  vertical  mixing  during  the  tests  was  negligible. 
The  rings  were  visible  continuously  during  the  linearly  strati- 
fied tests,  but  the  highly  refractive  interface  between  fresh 
water  and  brine  in  the  density  discontinuity  tests  obscured  the 
ring  during  its  passage  through  the  interface,  preventing  measure 
ments  of  ring  geometry  during  this  portion  of  the  test. 

The  vortex-ring  vertical  descent  z , diameter  D , and  the 
major  and  minor  semi-axes  (a  and  b,  respectively)  of  the  oval  of 
fluid  carried  along  by  the  ring  during  its  descent  were  measured 
from  projections  of  the  movie  film.  It  was  not  possible  to  ob- 
serve the  core  structure  or  velocity  with  this  method.  The  verti 
cal  descent  and  ring  diameter  were  relatively  easy  to  measure 
during  the  initial  portion  of  each  run,  but  measurement  of  ring 
diameter  became  progressively  more  difficult  during  the  latter 
portion  of  descent.  During  approximately  the  last  half  of  each 
run,  the  diameter  could  not  be  measured  due  to  turbulent  detrain- 

ment  of  dye.  The  oval  size  was  difficult  to  measure.  The  oval 
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was  measured  by  assuming  that  the  dye  carried  along  with  the  ring 
defined  the  oval  boundaries.  This  was  probably  not  true  in  the 
initial  ^/ortion  of  each  run,  but  diffusion  of  dye  during  descent 
appeared  to  fill  the  oval  sufficiently  to  provide  a meaningful 
measurement,  at  least  in  the  middle  portion  of  each  test.  In  the 
tests  carried  out  with  the  sharp  density  discontinuity,  refraction 
of  light  at  the  cell  boimdary  after  ring  penetration  of  the  heav- 
ier density  fluid  provided  a very  good  oval  definition  in  these 
cases  for  some  time  after  penetration  although  the  minor  axis  b 
was  still  not  as  well  defined  as  the  major  axis. 

The  Froude  number*  defined  for  the  density  discontinuity 
and  for  the  linear  gradient  cases  (and  given  in  Table  4-1)  differ. 
For  the  discontinuity 


where 


Ap/pQ  = 0.19 


^’^l  = 


and  R_ 


“gAp  ~ 

fo\ 


T7T 


(4.1) 


reaching  the  discontinuity. 


Fro  = 


is  the  initial  ring  radius  before 
For  the  linear  density  gradient 

^o 

inr  , (4.2) 


where 


N = - 


g dfTI' 


(4.3) 


is  the  Brunt-Vaisala  frequency,  equal  to  1.41  rad/sec  for  the 
present  tests. 


Shown  in  Table  4-1  is  the  maximum  depth  of  penetration  of 
the  rings,  (z^^^-z^) , measured  to  the  center  of  the  cell 
an  arbitrarily  selected  initial  point  during  descent. 


IS 

o 


Since  three  definitions  of  Froude  number  have,  through  necessity, 
been  used  in  this  report,  they  will  be  distinguished  by  the 
subscripts  1,  2,  and  3. 
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The  cell  size  for  the  non-stratif ied  cases  did  not  change 
appreciably  during  descent,  and  the  ring  reached  the  bottom  of 
the  test  tank.  In  the  stratified  cases,  the  oval  size  decreased 
appreciably,  and  in  all  stratified  cases,  the  rings  dispersed 
before  reaching  the  bottom  of  the  tank. 

Of  particular  note  was  the  structure  of  the  ring  during 
descent,  which  initally  appeared  very  ordered,  but  which  even- 
tually broke  down  due  to  an  instability  which  has  been  observed 
many  times  in  experimental  investigations  of  this  type.  See,  for 
instance,  Widnall,  Ref.  4,  for  a discussion  of  this  phenomenon. 

The  swirl  around  the  core,  which  could  be  seen  quite  clearly  during 
the  initial  stages  by  observing  the  entrained  dye  traces,  became 
impossible  to  see  at  the  latter  stages  of  descent  although  the  core 
region  was  still  visible.  Finally,  the  structure  was  completely 
obscured  by  refraction  produced  by  the  highly  turbulent  flow  of 
mixed  fluid  throughout  the  oval.  No  core  region  could  be  seen  at 
that  point,  and  the  ring  finally  was  destroyed  in  a catastrophic 
process  as  described  by  Maxworthy,  Ref.  16. 

4.2  COMPARISON  OF  DATA  AND  NUMERICAL  RESULTS 

The  comparison  of  experimental  data  with  the  numerical  re- 
sults of  the  stratified  axisymmetric  code  are  included  in  Figures 
4-1  through  4-13.  The  non-stratif ied  results  are  also  shown  in 
Figures  4-1  and  4-2.  Both  theory  and  experiment  reveal  no  dramatic 
behavior  for  the  non«stratified  case  and  the  ring  is  seen  to  propa- 
gate at  nearly  constant  velocity.  Figure  4-1,  while  increasing  in 
radius  only  slightly.  Figure  4-2.  The  horizontal  oval  axis  is 
predicted  well,  but  during  the  initial  portion  of  the  run,  the 
simulation  appears  poor  for  the  vertical  extent  of  the  oval.  This 
is  probably  because  the  dye  has  not  yet  had  time  to  diffuse  through 
the  oval  as  mentioned  previously.  Note  that  the  values  of  a/R^ 
and  b/R  obtained  from  the  constant  core  vorticity  assumption  are 

plotted  as  a dashed  line  in  Figure  4.2  for  the  vorticity  spread 
used  in  the  run . 

Figures  4-3a  and  4-3b  present  the  calculated  stream  function 
for  the  non-stratif ied  case  for  non-dimensional  times  of  2.0  and 
10.0.  These  plots  illustrate  that  little  if  any  change  is  pre- 
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dieted  in  the  oval  volume  for  the  non-etratif ied  case.  This  type 
of  plot  will  be  presented  later  for  the  stratified  cases,  where 
a very  noticeable  decrease  in  oval  size  will  be  shown.  Note  that 
the  abscissa  in  these  figures  indicates  the  vertical  position 
measured  from  the  initial  position  z^  . 

The  measured  and  calculated  trajectory  for  the  case  3 densi- 
ty discontinuity  ring  injection  is  given  in  Figure  4-4,  and 
the  radius  and  oval  parameter  histories  are  shown  in  Figure  4-5. 
Again  shown  in  the  figures  are  the  values  of  s/Rq  3nd  b/R^  for 

the  value  of  r /R  used  for  the  calculation.  In  this  case. the 
c o ’ 

measured  initial  value  of  ^/Rq  was  misleading  since  a vorticity 
spread  based  on  it  gave  a poor  simulation  and  the  case  was  rerun 
using  the  value  given  in  Table  4-1.  The  discrepancy  may  again 
be  due  to  the  failure  of  the  dye  to  define  adequately  the  oval 
boundary . 

Figures  4-6  and  4-7  are  the  calculated  stream  function  and 
density  isopleths  for  case  3.  At  tV^/R^  = 2.0,  the  ring  is 
approaching  the  density  discontinuity,  while  at  tV^/R^  = 6 , it 
has  fully  penetrated.  The  decrease  in  oval  size  in  this  case  is 
quite  evident.  Note  the  intensity  of  the  density  gradient  be- 
tween the  fluid  within  the  oval  and  the  ambient  fluid.  It  is 
also  interesting  that  there  is  little  predicted  difference  be- 
tween the  density  levels  and  the  shape  of  their  isopleths  in  the 
vicinity  of  the  core.  The  suppression  of  turbulence  in  the  core 
region,  as  a consequence  of  the  high  stability  of  the  vortex  flow, 
prevents  mixing  of  the  core  fluid  with  entrained  ambient  fluid 
in  the  outer  region  of  the  oval. 

The  linear  stratification  comparisons  are  given  in  Figures 
4-8  through  4-13.  The  ring  trajectories  for  cases  4,  5,  and  6 
are  shown  in  Figures  4-8,  4-9,  and  4-10,  respectively.  The  numer- 
ical simulations  for  these  cases  are  excellent  over  most  of  the 
descent,  with  the  possible  exception  of  the  lowest  injection 
velocity.  Figure  4-8.  It  was  noted  during  the  tests  that,  for 
injection  velocities  even  lower  than  those  analyzed,  the  rings 
were  not  completely  formed  at  initial  penetration  into  the  strati- 
fication and  tended  to  tilt  progressively  from  the  vertical. 
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The  prediction  of  ring  radius  and  oval  parameter  history  also 
gives  excellent  agreement  over  the  measurable  portion  of  descent, 
Figures  4-11  through  4-13.  The  poor  agreement  for  b/R^  is 
again  evident  in  the  initial  portion  of  the  run. 

Figures  4-14a  through  4-14d  present  the  calculated  stream 
function  for  the  linear  stratification  with  Froude  number  6.57, 
case  6.  Here,  the  decrease  in  oval  size  is  very  striking  at  the 
point  where  the  descent  has  almost  stopped. 

Figures  4-15  and  4-16  are  plots  of  the  ring  radii  versus  time. 
Figure  4-15  gives  the  measured  non-dimensional  radius  for  cases  2 
and  3,  which  are  ring  injections  into  a density  discontinuity. 

Over  the  range  of  time  in  which  the  radius  could  be  measured  after 
passing  the  discontinuity,  the  radius  history  exhibits  a t”^  de- 
pendence. Figure  4-16  compares  the  numerical  results  of  the  WAKE 
code  for  both  the  discontinuity  and  linear  stratification.  The 
calculated  result  for  the  discontinuity  is  very  close  to  a t ^ 
dependence.  For  a linear  stratification,  the  radius  decrease 
approaches  a t ^ dependence.  It  was  not  possible  to  verify  this 
behavior  during  the  experiment  since  the  radius  could  not  be 
measured  during  the  latter  stages  of  descent. 

Figure  4-17  is  a plot  of  the  non-dimensional  impluse 
P/PqR^  , kinetic  energy  circulation  I’/R^^V^ 

for  case  6,  a linear  stratification  run.  The  kinetic  energy  does 

not  appear  to  exhibit  any  particular  functional  dependence  on 

time,  but  over  the  middle  portion  of  the  calculation,  the  impulse 
-2  -S 

behaves  as  t and  the  circulation  as  t ^ . These  calculations 
are  based  on  integrations  over  the  entire  numerical  grid  however, 
and  therefore  the  plotted  results  include  the  contribution  of 
the  countersign  vorticity  generated  at  the  oval  boundary.  It  is 
noted  that  Saffman,  Ref.  7,  has  shown  that  the  concept  of  impulse 
in  a stratified  fluid  is  not  generally  useful  unless  this  counter- 
sign vorticity  is  included  and  the  Boussinesq  approximation  is 
made.  Both  criteria  are  met  here. 
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The  objective  of  the  experimental  and  numerical  simulation 
comparison  reported  in  this  section  was  to  provide  verification 
of  the  A.R.A.P.  second-order  closure  turbulent  model.  The 
agreement  between  the  numerical  results  based  on  the  model  and 
the  experimental  data  has  provided  this  verification. 
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Figure  4-1.  Vortex-ring  trajectory  for  non-stratif ied 
fluid,  case  1.  a Data;  — axisymmetric 
wake  code  calculation. 
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Figure  4-3b.  Stream  function  for  ease  1 at  t V K » 10. 0 
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Figure  4-4.  Vortex-ring  trajectory  for  penetration  of  density 
discontinuity,  case  3.  o Data;  — axisymmetric  stratified 
WAKE  code  calculation. 
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»-6a.  Stream  function  for  vortex  ring  descending  into  a 


(z-c  = 2.0  at 
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density  discontinuity,  case  3. 

2.0.  Stream  function  contours  from 
-0.6  to  d.*!  at  intervals  of  0.3. 
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Figure  4-6b.  Stream  function  for  case  3.  (z-z^)/Rq  = 5.4  at 

tV  /R  = 6.0  . Stream  function  contours  from 
o o 

^ -0.4  to  +3.2  at  intervals  of  0.2. 
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Figure  4-7a.  Isopyenals  for  vortox-ring  ile«cent  into  a dt'nslty  '' 

discontinuity,  case  3.  Ring  starling  penetration  . 

of  discontinuity,  tV  /R  - 2.0,  (z-z  Wr  ■"  2.0.  f 

’ o o • ' o o 

Contours  front  -0.4b  to  -H).4b  at  Inlerv.tls  itf  0.1b. 
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MAXIMUM  MEASURED  PENETRATION 


calculation 
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figure  4-10.  Vortex-rlng  trajectory  In  a linearly  atratifled  fluid  for  Fr^  - 6.57  (caae  6) 
a Data;  — axisynmetric  stratified  WAKE  code  calculation. 
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Figure  4-12.  Vortex-ring  semi-major  axis  history  for  Fr.,  “ I>.18  , 

b.l7  and  6,57  teases  4,  5 and  6).  o Data;  — axisymmotr ic 
stratified  WAKF  code  ealeulntion;  — — — non-st rat  if led 
constant  core  vorticity  value. 
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5.  COMPARISON  OF  IW- DIMENSIONAL  WAKE  CALCULATIONS  WITH  FULL-SCALE 

FLIGHT  DATA 

5.1  TWO-DIMENSIONAL  WAKE  CODE 

The  code  validation  contained  in  our  vortex  ring  study  pro- 
vides increased  confidence  in  the  two-dimensional  results  to  be 
compared  with  flight  data.  For  vortex  flows  having  moderate  vor- 
ticity  gradients,  the  solution  procedure  of  the  two-dimensional 
equations,  which  are  included  in  Appendix  A,  may  be  identical  to 
that  used  in  solving  the  axisymmetric  equations.  However,  due  to 
the  severity  of  the  vortex  wake  chosen  for  the  calculations  to  be 
presented  here,  a new  technique  for  calculating  strong  vortex  wakes 
has  been  developed.  Details  of  this  technique  will  be  presented 
in  following  sections. 

5.2  DATA  SELECTION 

To  reduce  extraneous  effects  which  cannot  be  included  in  a 
numerical  calculation,  criteria  were  set  up  as  a guide  in  selecting 
the  flight  case  to  be  studied.  A measurement  made  on  a clean 
aircraft  configuration  (minimum  flaps,  no  spoilers)  in  level  flight  j( 

was  sought,  and  calm  ambient  conditions,  e.g.,  low  cross  shear 
and  ambient  turbulence,  were  factors  considered.  The  iiumber  and 
quality  of  measurements  made  during  the  test  were  also  major  factors 
in  selection.  Finally,  of  course,  stable  atmospheric  stratification 
was  a required  condition. 

The  measurement  selected  for  comparison  with  our  numerical 
simulation  was  Run  8 of  the  1975  Rosamond  Dry  Lake  tests,  made 
using  a Boeing  747  aircraft.  A complete  description  of  these  tests 
and  data  for  all  runs  have  been  reported  by  Burnham,  et  al..  Ref.  30. 

The  measurements  of  vortex  characteristics  were  taken  with  tlie  air- 
craft in  a variety  of  flap  settings,  landing  gear  and  spoiler  f 

configurations,  and  meteorological  conditions.  Takeoff,  level 
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flight,  ami  landing  aircraft  attitudes  also  were  varied.  A 
total  of  54  passes  were  flown  at  altitudes  varying  from 
30-  to  244-  ni.  Test  instrumentation  included  a nionostatic 
acoustic  system  for  vortex  tracking  and  strength  determination, 
three  anemometer  baselines  for  vortex  tracking,  and  a mobile  Laser- 
Doppler  Veiocimeter  (LDV)  for  collecting  detailed  velocity 
signatures . 

5.3  ThST  CONDITIONS 

The  important  flight  and  environmental  pariimeters  for  Run  8 

of  the  1975  Rosamond  Dry  Lake  tests  is  given  in  Table  5-1,  taken 

from  Table  2-3  of  Ref.  30.  The  total  circulation  of  the  B-747  air- 

2 

craft  for  this  test  run  was  reported  to  be  628  m^/sec,  the  initial 
vortex  separation  to  be  50  m , and  the  initial  descent  rate  to 
be  2 m/sec.  A total  of  about  80  sec  of  time  test  was  recorded 
for  this  run. 

The  measured  swirl  velocity  profile  for  Run  8 at  11  seconds 
alter  passage  of  the  generating  aircraft  is  shown  in  Figure  5-1, 
and  tlie  associated  circulation  profile  is  given  in  Figure  5-2. 

The  curves  through  the  data  are  based  on  a constant  vorticity  core 
and  a Hof fman-Joubert  turbulent  vortex  circulation  profile  (Ref.  31), 
beyond  a radius  of  1.5  m . At  10.0-m  radius,  the  velocity 
is  assumed  to  depend  upon  the  inverse  radius  (zero  vorticity). 

The  swirl  velocity  as  shown  in  Figure  5-1  represents  a severe 
vortex  wake,  with  a maximum  swirl  velocity  of  nearly  17  m/sec 
occurring  at  a radius  of  2.5  m . To  resolve  this  profile  numeric- 
ally requires  mesh  points  in  the  computation  at  less  than  1-m 
intervals.  Because  we  must  also  resolve  the  effects  of  stratifi- 
cation, which  are  occurring  at  the  edges  of  the  oval  more  than 
25  m from  the  core  region,  and  to  provide  a large  enough  domain 
to  avoid  wave  interference  effects  at  the  computational  boundaries, 
the  outer  boundary  (horizontally)  must  be  placed  at  least  at  75  m 
from  the  wake  centerline,  wltich  because  of  symmetry  constitutes  one 
of  the  vertical  computational  boundaries.  The  WARK  code  vs  pres- 
ently limited  to  a 33x40  mesh  on  our  DKC®  l’DF-11/70  computer,  and 
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TABLE  5-1.  FLIGHT  PARAMETERS  AND  METEOROLOGICAL  FACTORS 
OF  RUN  8,  1975  ROSAMOND  DRY  LAKE  TESTS  (REF.  30) 


! 


Flight  time 

0745  PST 

Altitude 

244  m 

Airspeed 

77  m/sec 

Aircraft  weight 

244000  kg 

Flap  setting 

30/30  (deg) 

Landing  gear 

Down 

Spoilers 

None 

Thrust  (EPR) 

1.20 

Wind  speed 

1 m/sec 

Wind  direction 

310  deg 

1/3  * 

Atmospheric  dissipation  rate  (e  ) 

1.4  cm  /sec 

Lapse  rate 

1 °C/100  m 

Flight  attitude 

Level 

*At  flight  altitude. 


therefore,  a calculation  requiring  the  resolution  of  a small-scale 
phenomenon  on  a large  numerical  domain  is  beyond  the  capabilities 
of  the  code,  even  if  a nonuniform  grid  is  employed.  If  the  code 
were  reconfigured  to  allow  more  grid  points  (necessitating  a far 
greater  amount  of  disk  I/O),  a time  step  of  0.05  seconds  would 
be  needed  in  order  to  satisfy  the  CFL  (CouranL-Friedr ichs-Lewy ) 
condition  due  to  the  large  swirl  velocity.  To  simulate  70  sec  of 
test  time  would  require  approximately  1400  time  steps,  which  would 
use  more  than  100  hours  of  CPU  (central  processor  unit)  time  on 
our  PDP-11/ 70, and  at  least  three  times  as  much  real  time  due  to 
the  excessive  number  of  disk  operations  required  for  such  a fine 
grid. 

The  feasibility  of  the  computation  would  be  in  doubt  if  the 
calculation  were  carried  out  along  traditional  lines.  A new 
approach  was  therefore  initiated  for  the  computation  as  described 
below. 

5.4  VORTEX  WAKE  CALCULATION  BASED  ON  A SPLITTING  TECHNIQUE 

From  an  examination  of  the  data,  it  was  noticed  that  the  de- 
cay in  the  circulation  profile  was  small  over  most  of  the  test 
period.  It  was  decided  that  the  departure  of  the  vortex  from  its 
initial  state  was  the  basis  of  an  improved  calculation  procedure. 
The  computation  would  consequently  be  relieved  of  supporting  the 
advection  of  vorticity  and  turbulence  by  a purely  numerical  cal- 
culation of  the  advection  velocities. 

If  the  vorticity  were  written  as  the  sum  of  the  initial  vor- 
ticity distribution  and  a time-varying  perturbation  about  the 
initial  state  (hence  the  designation,  splitting  technique),  and 
substituted  into  the  governing  equations,  we  would  obtain  equations 
which  may  be  solved  for  the  time-dqpendent  portion.  Perturbation 
about  the  initial  state  does  not  imply  an  approximation  here  in 
tlie  sense  of  small  variations.  The  procedure  is  an  exact  one  and 
is  done  to  reduce  the  magnitude  of  the  numerically  calculated 
vorticity  and  velocities.  Depending  upon  the  complexity  of  the 
analytical  expression  selected  Cor  the  initial  vorticity  profile, 
the  initial  conditions  can  be  represented  as  closely  as  we  choose. 


The  benefits  gained  by  this  procedure  are  twofold.  First,  the 
analytical  vorticity  and  velocities,  which  will  be  large  compared 
to  perturbation  quantities  over  most  of  the  calculation  time,  are 
invariant  with  time  at  each  grid  point.  Numerical  noise  is  there- 
fore reduced  considerably.  Second,  the  perturbations  are  initially 
zero,  and  if  they  are  to  remain  small  over  the  computation  life- 
time, the  perturbation  gradients  will  also  be  small.  This  relaxes 
somewhat  the  mesh  size  requirements,  opens  up  the  allowable  time 
step,  and  simultaneously  reduces  the  risk  of  numerical  instability. 

To  provide  a smooth  generating  function  for  the  initial  vor- 
ticity, a Gaussian  profile  is  chosen  to  match  the  swirl  velocity 
and  circulation  profiles  as  closely  as  possible  with  this  type  of 
distribution.  The  vorticity  used  for  the  calculations  is 

r-(r,/3.7)2  -(r.,/3.7)2"l 

= 14.4  l^e  ^ - e " J (5.1) 

in  sec”^,  where  r^^  is  measured  from  y - 25,  z = 0 m,  and  V2 
from  y - -25,  z = 0 m.  The  swirl  velocity  and  circulation  pro- 
file derived  from  this  vorticity  are  compared  with  data  in 
Figures  5-3  and  5-4.  These  profiles  do  not  represent  the  measured 
velocity  and  circulation  as  closely  as  the  profiles  of  Figures  5-1 
and  5-2.  They  are  generated  by  three  different  vorticity  functions 
over  the  radius  of  the  vortex,  and  introduce  discontinuous 
derivatives  at  two  regions  of  the  calculation.  The  Gaussian  pro- 
file is  used  to  avoid  this  possible  source  of  numerical  instability. 
Note  that,  since  we  are  calculating  the  perturbation  vorticity,  it 
should  be  possible  to  use  a generating  function  which  is  only 
approximately  equal  to  the  actual  profile.  If  our  calculation 
does  indeed  simulate  the  physics  of  the  vortex  flow,  the  generating 
function  plus  perturbation  should  evolve  toward  the  measured  vortex 
profiles  as  the  computation  proceeds. 

To  accommodate  the  perturbation  procedure,  the  mean  vorticity 
equation  in  Cartesian  coordinates  may  be  written  in  the  following 
form ; 
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where  the  zero  subscript  here  denotes  the  analytically  derived  ve- 
locities due  to  the  initial  vorticity  field.  The  perturbation 
vorticity  and  velocities  are  denoted  by  primes.  The  perturbation 
stream  function  is  calculated  from 
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The  vortices  were  kept  in  the  center  of  the  computational 
domain  by  an  algorithm  which  traced  the  maximum  total  vorticity 
t ” ^ kept  this  at  z = 0 , by  adding  an  appropriate 

vertical  washing  velocity  into  the  calculation. 

The  perturbation  and  initial  velocities  entered  into  the  calcu- 
lation of  the  turbulence  correlations  in  both  the  advection  terms 
and  the  velocity  gradient  terms.  The  advection  terms  in  the 
Reynolds  stress  equations  were  treated  as  shown  above  for  the 
vorticity  equation.  The  velocity  gradients  in  the  production 
terms  of  the  turbulence  equations  were  written  as  the  sum  of  the 
analytically  derived  gradients  of  the  initial  velocity  and  the 
numerically  calculated  gradients  of  the  perturbation  velocities. 

A shortwave  filter  algorithm  (Sheng,  Ref.  32)  was  used  to  control 
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numerical  noise  on  the  perturbation  vorticity  and  turbulence  velocity 
correlations.  However,  the  filter  was  used  only  when  the  criterion 
on  maximum  shortwave  oscillations  (local  peaks  and  curvature)  was 
exceeded  in  two  directions  simultaneously. 


5 . 5 AMBIENT  STRATIFICATION 


The  stable  ambient  stratification  at  the  flight  altitude  is 
measured  to  be  1 °C/100  m.  During  descent,  a vortex  pair  compresses, 
however,  at  a rate  of  approximately  1 °C/100  m if  it  is  assumed 
the  compression  is  adiabatic.  The  calculations  are  performed  using 
potential  temperature  as  the  independent  variable  and  the  stratifi- 
cation set  equal  to  2 °C/100  m.  The  effective  Brunt-Vaisala 
frequency  is  therefore 


Gt)3i 


= 0.0265  rad/sec 


where  0 = 280  °C  . Defining  the  Froude  number  to  be 


Fra  = 


r 

o 



2Trs^N 


(5.5) 


(5.6) 


the  Froude  number  for  Run  8 at  11  seconds  after  aircraft  passage  is 
equal  to  6.03. 


5.6  INITIAL  AND  AMBIENT  TURBULENCE 

The  initial  turbulence  within  the  vortices  could  not  be 
measured.  An  estimate  of  the  initial  turbulence  was  made  by  em- 
ploying the  observed  time  required  for  smoke  entrained  in  the  vortex 

to  diffuse  to  the  edge  of  the  cell,  as  reported  in  Ref.  30.  This 

2 2 

led  to  an  initial  turbulence  level  q of  about  1 m /sec 

The  ambient  turbulence  may  be  derived  from  the  measured 
ambient  dissipation  rate  e . Based  on  an  assumed  ambient  turbu- 
lence scale  of  2.5  m , the  ambient  turbulence  is  about 
2 2 

0.03  m /sec  . This  level  should  not  appreciably  affect  the 

behavior  of  the  vortices  over  the  time  scale  of  the  self-generated 
turbulence  and  stratification  effects  on  the  vortices. 

65 


5.7  TURBULENCE  SCALE 


The  macroscale  A of  the  self-generated  turbulence  is  one 
of  the  most  uncertain  variables  of  the  vortex  flow  phenomenon  we 
are  investigating.  It  is  expected  that  this  scale  may  be  controlled 
by  the  characteristic  dimensions  of  the  flow  field,  which  in  this 
case  are  three  in  nxxmber;  the  core  radius,  the  core  separation, 
and  the  length  associated  with  the  gradient  of  density  occurring 
at  the  boundary  of  the  cell  of  recirculating  fluid.  Generation 
of  the  turbulence  and  the  length  scale  of  this  turbulence  is  a 
function  of  the  core  size  and  the  density  gradient  scale  in  the 
respective  regions  of  turbulence  production.  Using  the  vortex-ring 
results  as  a guide,  where  “ 0.03  gives  the  best  results. 

A/s  • 0.03  is  used  in  the  present  analysis. 

5.8  COMPARISON  OF  WAKE  CALCULATIONS  WITH  DATA 

The  data  reported  by  Burnham  et  al.  in  Ref.  30  included 

histories  of  the  vortex  pair  descent  velocity  and  altitude, 

separation  between  vortices,  maximum  tangential  velocity,  and  the 

vortex  circulation  and  average  circulation  vs.  radius.  The  average 

circulation  vs.  radius  was  defined  in  Ref.  30  as 

r 

r(r)dr  (5.7) 

The  calculated  vortex-pair  descent  velocity  and  altitude  are 
compared  with  the  measured  values  in  Figures  5-5a  and  5-5b,  re- 
spectively. The  agreement  between  measurement  and  calculation  is 
excellent.  Note  that  the  measured  descent  velocity  and  vortex 
altitude  oscillate  about  a mean  which  is  represented  by  the 
calculated  values.  These  oscillations  are  similar  to  those  that 
one  would  expect  to  see  by  observing  the  spiriling  motion  of  one 
vortex  about  another.  This  phenomenon  did  occur  during  the  tests 
where  inboard  flap  generated  vortices  were  seen  to  spiral  about 
and  merge  with  the  tip  vortices.  The  period  of  rotation  for 
this  interaction  was,  however,  considerably  shorter  than  30  sec- 
onds. and  merging  of  these  vortices  occurred  about  3 wing  spans 
behind  the  aircraft. 
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Note  that  the  calculated  tip  vortex  descent  velocity  decreases 
monotonically  and  does  not  accelerate  as  predicted  by  several 
previous  analyses,  see  Table  1-1. 

Figure  5-6  shows  a comparison  between  the  measured  and  calcu- 
lated average  circulation,  defined  by  Eq.(5.7),  at  5-,  10- , and 
15- , m radii  from  the  center  of  the  core.  ^ avg  calculated 

using  a radially  varying  r(r)  , which  in  turn  was  calculated 
based  on  an  area  integration  of  the  vorticity  out  to  radius  r . 

The  point  r = 0 corresponded  to  the  maximum  stream  function. 

The  integration  was  carried  out  in  a circle  with  maximum  radius 
defined  by  the  distance  from  r = 0 to  the  centerline  between  the 
vortices . 

The  profiles  of  circulation,  r(r)  vs . r , from  11  to  76 
sec  are  presented  in  Figures  5-7a  to  5-7f.  The  characteristics 
of  the  measured  circulation  and  average  circulation  profiles  are 
represented  by  the  calculations.  The  initial  rate  of  decrease 
in  r and  'ayg  small,  increasing  gradually  to  the  end  of  the 
measurement.  The  decrease  in  circulation  with  radius  as  the  edge 
of  the  recirculation  cell  is  approached  is  due  to  the  countersign 
vorticity  in  this  region.  This  decrease  is  predicted  by  the 
calculation,  and  can  be  seen  in  the  data. 

Figure  5-8  shows  a comparison  between  the  measured  maximum 
tangential  velocity  and  the  calculated  value.  Again,  the  agree- 
ment is  excellent  over  most  of  the  data  period.  Note  that,  even 
after  the  vortex  pair  has  stopped  descending,  there  is  still  an 
appreciable  swirl  velocity  in  the  core.  The  vortices,  therefore, 
have  considerable  strength  remaining  although  they  have  been 
stopped  in  their  descent  by  the  stratification.  This  result  has 
been  previously  found  in  our  atmospheric-wake  vortex  interaction 
study  described  in  Ref.  13. 

The  separation  between  the  vortices  is  shov.ni  in  Figure  5-9. 

The  data  here  extend  out  only  to  about  25  sec,  when  one  of  the 
vortices  drifts  beyond  the  scanner  elevation  angle  limit . The 
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separation  is  calculated  to  remain  nearly  unchanged  over  more  than 
one-half  of  the  descent  time,  decreasing  by  more  than  10  percent 
of  the  initial  separation  only  during  the  final  stages  of  descent. 
This  decrease  in  separation  is  more  than  found  previously,  see 
Bilanin  et  al . , Ref.  13.  The  vortex  studied  there  is  an  initially 
more  diffuse,  aged  vortex,  while  the  present  one  is  in  a much 
earlier  stage.  The  aged  vortex  is  dominated  by  turbulent  diffusion 
in  and  about  the  core  region.  Turbulence  production  in  the 
center  of  the  core  of  the  measured  vortex  is  suppressed  by  the  cen- 
trifugal effects  of  the  swirling  flow,  and  turbulence  there  is  due 
only  to  that  diffused  inward  or  due  to  axial  gradients  of  velocity, 
buring  the  time  in  which  stratification  effects  are  occurring,  the 
present  vortices  are  dominated  by  the  turbulence  in  the  cell  bound- 
ary region.  Figures  5-lOa  to  5-lOd  are  contour  plots  of  the  vor- 
ticity  in  the  core  region,  and  show  a very  slow  decay  over  most  of 
the  test  period.  Note  that  the  vertical  coordinate  z in  these 
and  following  figures  is  measured  from  the  vortex  pair  altitude  at 
11  sec  after  passage  of  the  aircraft,  220  m. 

Figure  5-11  shows  contours  of  the  calculated  vorticity  at  65 
sec  in  the  core  region,  with  contour  intervals  of  3 sec"^,  and  in 
the  diffuse  region  containing  the  countersign,  stratification 
generated  vorticity,  with  intervals  of  0.05  sec"^.  Although  small 
in  value  locally,  the  large  extent  of  the  countersign  vorticity  is 
sufficient  to  stop  the  descent  of  the  primary  vorticity  in  and 
around  the  cores.  A calculation  of  the  descent  velocity  based  on 
the  strength  and  centroid  of  the  primary  and  countersign  vorticity 
at  65  sec  predicts  a descent  velocity  of  0.9  m/sec.  The  actual 
mean  descent  velocity  of  the  numerical  calculation  is  0.4  m/sec. 

This  disparity  is  due  to  the  very  diffuse  extent  of  the  counter- 
sign vorticity,  and  a simple,  concentrated  vorticity  argument 
cannot  be  used  accurately.  The  simplified  calculation  does,  however, 
predict  a decrease  in  vortex-descent  velocity  and  not  an  increase  as 
occurs  if  only  the  core  separation  is  used  as  a basis  for  the  cal- 
culation . 
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The  balance  between  the  primary  and  countersign  vorticity  can 
be  appreciated  by  seeing  their  variation  with  time,  Figure  5-12, 
which  is  plotted  non-dimensionally . The  total  circulation  r 
decreases  by  nearly  one-half,  a consequence  of  the  countersign  cir- 
culation . The  circulation  of  the  primary  vortex,  = r - 1’^ 

remains  high  however.  Note  that  is  nearly  linear  in  time  for 

this  vortex.  Also  shown  in  Figure  5-12  are  T , the  total  kinetic 
energy  of  both  vortices,  which  is  initially  equal  to  that  of  an 
elliptically  loaded  wing,  and  P , the  impulse.  There  appears  to 
be  no  simple  time  dependence  for  either  the  energy  or  the  impluse. 

The  calculated  streamline  patterns  at  11,  39,  55,  and  65  sec 
are  shown  in  Figures  5-13a  through  5-13d,  respectively.  The 
recirculation  cell  area  decreases  gradually  toward  the  end  of 
descent,  again  in  agreement  with  our  previous  findings  reported 
in  Ref.  13. 

2 

The  last  figures,  5-14a  to  5-14c,  are  contour  plots  of  q 
at  18,  39,  and  55  sec,  respectively,  showing  the  development  of 
the  turbulence  field  with  time.  The  level  of  turbulence  is  a 
maximum  in  a ring  of  approximately  6-m  radius,  decreasing  to 
10  percent  of  this  maximum  at  the  center  of  the  core.  After  the 
field  has  developed,  it  remains  remarkably  constant  with  time  due 
to  the  very  tight,  stable  core  structure.  The  maximum  velocity 
associated  with  the  turbulence  is  approximately  10  percent  of  the 
maximum  swirl  velocity. 

No  temperature  measurements  were  made  for  the  tests  analyzed 
here.  Most  notable  about  the  calculated  temperature  field  was  the 
uniformity  of  temperature  throughout  the  recirculation  cell.  The  po 
tential  temperature  of  the  cell  remained  very  nearly  at  its  initial 
value,  indicating  only  very  slight  entrainment  of  ambient  fluid 
during  descent. 


Tangential  velocity,  m/sec 


r 


Figure  5-1.  Measured  Boeing  7A7  vortex  swirl  velocity  vs.  radius 
at  11  sec  after  aircraft  passage  and  analytical 
representation  based  on  constant  vorticity  core 
(r  < 1.5m),  Hof fman-Joubert  turbulent  profile 
(1.5  <r  < 10m),  and  irrotational  velocity  profile 
(r  s lotn). 
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FiKurc  5-4.  Measured  circulation  distribution  and  analytical 
profile  based  on  a Gaussian  vorticity  distribution 
at  11  sec. 
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Figure  5-5.  Comparison  between  measured  and  calculated  vortex 
descent  velocity  (a)  and  altitude  (b)  vs.  time. 
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-7a.  Comparison  between  measured  and  calculated  initial 
circulation  at  11  sec  after  aircraft  passage. 
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Figure  5-7c.  Comparison  between  measured  and  calculated 

circulation  at  25  sec  after  aircraft  passage. 
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Figure  5-7d.  Comparison  between  measured  and  calculated 

circulation  at  39  sec  after  aircraft  passage 
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5-7f.  Comparison  between  measured  and  calculated 

circulation  at  76  sec  after  aircraft  passage 
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Figure  5-8.  Comaprison  between  measured  and  calculated  maximum 
swirl  velocity  vs.  time. 
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Figure  5-9.  Measured  and  calculated  vortex-pair  separation  vs. 
time . 
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Figure  5-lOa.  Contours  of  vorticity  in  the  core  region  at  11  sec. 
Contour  interval  2 sec”^. 
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Figure  5- 10c.  Contours  of  vorttc 
Contour  interval  2 
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Figure  5“L3a.  Streamline  pattern  at  11  sec.  Contour  interval 
of  20  m2/sec.  Descent  velocity  - 2.0  m/sec. 
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Figure  5-13d.  Streamline  pattern  at  65  sec.  Contour  interval 
of  20  m2/sec.  Descent  velocity  =0.4  m/sec. 
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Figure  5-14b.  Contours  of  turbulence  at  39  sec.  Contour  interval 

1 m2/sec2. 
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Figure  5-14c.  Contours  of  turbulence  at  55  sec.  Contour  interval 

1 m^/sec2. 


6.  DISCUSSION  AND  CONCLUSIONS 

6.1  DISCUSSION 

In  this  study,  we  have  closely  simulated  the  behavior  of  a 
turbulent  vortex  descending  in  a stable  density  stratification, 
showing  agreement  between  the  calculated  and  observed  behavior. 

What  we  find  is  a very  strong,  tight  vortex  which  decays 
slowly  during  descent,  and  which  has  appreciable  strength  after 
it  has  stopped  descending.  This  investigation  suggests  that  the 
core  vorticity  structure  remains  well  ordered  throughout  most  of 
the  descent,  while  the  vorticity  within  the  core  is  slowly  diffused 
outward  due  to  a ring  of  turbulence  surrounding  the  core.  The 
turbulence  reaches  a peak  at  a point  about  twice  the  radius  of 
maximum  tangential  velocity.  A portion  of  the  vorticity  diffused 
away  from  the  core  is  annihilated  by  countersign  vorticity  generated 
at  the  boundary  of  the  recirculation  cell,  part  of  which  is  itself 
diffused  inward  toward  the  core  by  the  turbulence  generated  at  the 
cell  boundary.  The  countersign  vorticity  remains  small,  less  than 
1 percent  of  the  magnitude  of  the  maximum  core  vorticity,  because 
it  is  being  continually  diffused  into  the  cell  or  swept  out  of 
the  vicinity  of  the  cell  by  descent.  Although  the  level  of  the 
countersign  vorticity  is  small,  the  area  of  the  wake  containing  it 
reaches  an  area  comparable  with  that  of  the  recirculation  cell. 

The  location  of  the  countersign  vorticity  is  primarily  above  and 
outboard  of  the  cores.  This  will  tend  to  drive  the  cores  together, 
which  is  seen  in  Figure  5-9,  but  the  net  effect  of  the  strength 
and  location  of  the  countersign  vorticity  with  respect  to  the  cores 
is  to  halt  descent.  The  circulation  of  the  countersign  vorticity 
reaches  a magnitude  of  nearly  50  percent  of  the  circulation  of 
the  initial  primary  vorticity,  while  the  primary  vorticity  is  re- 
duced by  only  about  20  percent  below  its  initial  value  at  the  end 
of  descent.  The  cores,  then,  remain  very  strong,  and  a vortex  pair 
cannot  be  assumed  to  have  dissipated  when  it  stops  descending. 


The  turbulence  near  the  center  of  the  core  is  repressed  due 
to  the  stability  of  the  flow  in  this  region,  and  remains  small  over 
the  descent  time  time  of  the  pair.  It  reaches  a level  of  only  10 
percent  of  the  peak  value  outside  the  core,  accounting  for  the 
persistence  of  the  vortices.  However,  the  rate  of  annihilation  of 
vorticity,  and  the  extent  and  strength  of  the  countersign  vorticity, 
are  direct  consequences  of  the  turbulent  diffusion  process. 

The  effect  of  stable  atmospheric  stratification,  then,  is 
found  to  halt  descent,  but  not  necessarily  to  completely  disorganize 
the  vortex  pair.  Increasing  the  stable  stratification  would  de- 
crease the  total  descent  distance,  although  no  quantitative  state- 
ment is  possible  here  since  we  have  analyzed  a single  stratifica- 
tion condition.  In  Ref.  13  the  descent  distances  of  vortices 
having  considerably  less  strength  than  analyzed  here  were  found  to 
be  directly  proportional  to  Froude  number.  This  may  not  be  indica- 
tive of  the  behavior  of  stronger  vortices,  however. 

6 . 2 CONCLUSIONS 

The  results  of  this  study  may  be  summarized  as  follows: 

a)  Turbulence  effects  must  be  accounted  for  to  arrive  at  an 
accurate  picture  of  vortex-pair  behavior  in  a stratified  flow. 

b)  Unalleviated  vortices  remain  strong  during  descent  due 
to  suppression  of  turbulence  within  their  cores. 

c)  The  vortices  stop  descending  due  to  a diffuse  region  of 
stratification-produced  countersign  vorticity  primarily  outboard 
and  above  the  cores . 

d)  The  core  separation  does  not  control  the  descent  velocity 
to  the  extent  predicted  by  several  previous  investigators. 
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APPENDIX  A 

SECOND- ORDER  CLOSURE  MODELED  EQUATIONS 

A.l.  EQUATIONS  IN  AXISYMMETRIC  COORDINATES 

The  modeled  axisymmetric  equations  of  motion  for  the  mean 
variables  and  turbulent  correlations  (denoted  by  < >)  may  be 
written  as  follows.  The  tilde  indicates  variables  that  have 
normalized  by  r to  insure  analyticity  at  r = 0. 
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(A. 18) 


u and  w are  found  by  solving  the  Poisson  equation  for  the 
stream  function  iJj 
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The  modeling  constants  are  assigned  the  values;  v = 0.3  ; 


A = 0.75  ; a - 2.5  ; b = 0.125  ; s = 1.8  ; s,  - -0.35  ; s. 


Sj  = 0.375  and  s^  = 0.8  . 


A. 2 EQUATIONS  IN  CARTES lAI^  COORDINATES 

In  Cartesian  coordinates,  the  equations  for  the  two- 
dimensional  WAKE  calculations  are  written  as  follows  : 
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APPENDIX  B 


VORTEX- RING  NUMERICAL  BOUNDARY  CONDITIONS 

In  Reference  23,  Bilanln  et  al.  devised  a new  technique  to 
specify  computational  boundary  conditions  when  vorticity  is  not 
extensively  distributed  throughout  an  infinite  two-dimensional 
fluid.  This  same  technique  has  been  employed  to  obtain  the  numer- 
ical boundary  condition  on  the  stream  function  for  the  incompress- 
ible axisymmetric  case.  Consider  the  axisymmetric  vorticity 
distribution  t of  Figure  B1  at  time  t . The  stream  function 
is  computed  from 


Z 


Figure  Bl.  Representation  of  vorticity  distribution. 
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The  solution  to  this  equation  may  be  written 
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where  the  denominator  of  the  integrand  is  the  length  s between 
a point  at  (r.z)  and  that  at  (r',z)  . and  1/(Atts)  is  the 
appropriate  Green's  function.  0 is  the  included  polar  angle 
between  these  points. 

The  radial  and  axial  velocities  are  obtained  from 
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Although  the  integral  defining  could  be  carried  out  to 
determine  the  stream  function  on  the  computational  boundary,  a 
computation  grid  containing  JxK  points  would  require  approximate- 
ly 3(J+K)/2  integrations  over  three-dimensional  space.  This 
would  be  prohibitive  in  terms  of  computer  time,  and  we  took  an 
alternative  approach. 

Selecting  our  numerical  computational  boundaries  so  that 
they  are  away  from  and  surround  the  region  containing  vor- 
ticity,  Eq . (B.2)  may  be  approximated  by  expanding  for  r >>  r' 
z >>  z'  , resulting  in 
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where 


r 


perturbs  the  density  field,  the  inflow  density  at  z_  is  not 
precisely  a linear  function  of  z . In  steady-flow  coordinates 
the  velocity  along  streamlines  is  reduced  below  the  self-induced 
ring  velocity,  V , and  approaches  V as  R2  -►  »>  . Referring 
to  Figure  B2, 


Figure  B2 . Effect  of  descending  vortex  ring  on  isopycnals. 


For  a linear  stratification  , 
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the  functions 


for  m and  n up  to  2 are  given  by 
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APPENDIX  C 

REPORT  OF  NEW  TECHNOLOGY 


The  work  under  this  contract  has  contributed  to  increased 
understanding  of  aircraft  wake  vortex  behavior  in  stable  stratifi- 
cation, as  well  as  to  providing  confidence  in  numerical  simulation 
techniques  in  the  prediction  of  turbulent  vortex  pairs  and  vortex 
rings  in  stratified  fluids. 
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